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This  report  was  prepared  by  Malcolm  Mellor,  Research  Physical.  Scientist, 
of  the  Experimental  Engineering  Division,  U.S.  Army  Cold  Regions  Research  and 
Engineering  Laboratory.  The  work  was  performed  as  part  of  DA  Project 
4A762730AT&2 ,  Design,  Construction,  and  Operations  Technology  for  Cold  Regions 
Task  Area  CS,  Combat  Support,  Work  Unit  029,  Explosives  and  Projectile  Impact 
Under  Winter  Conditions. 

The  author  thanks  R*  Gerard  and  P.  Sellmann  for  technical  review  of  the 
manuscript  - 
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DERIVATION  OF  GUIDELINES  FOR  BLASTING  FLOATING  ICE 
Malcolm  Melior 

Introduction 

In  1972,  guidelines  for  blasting  floating  ice  sheets  were  developed  by 
analyzing  all  the  test  data  that  were  then  available  (Melior  1972).  A 
decade  later,  the  1972  curves  were  revised  after  addition  of  new  data  to 
the  file  (Melior  1982)  .  Since  then,  a  few  more  data  sets  have  been  ob¬ 
tained,  and  the  curves  have  been  updated  again. 

In  developing  the  original  design  curves,  the  primary  concern  was 
practical  convenience,  and  certain  simplifying  assumptions  were  made.  It 
now  appears  that  these  assumptions  are  adequately  justified  within  broad 
limits  of  practical  applicability,  and  the  empirical  correlations  are  not 
improved  by  more  involved  treatment  of  the  relevant  physics. 

One  purpose  of  this  paper  is  to  give  the  current  version  of  the  ice¬ 
blasting  curves,  and  to  describe  their  use.  A  second  purpose  Is  to  explain 
alternative  methods  of  analysis  which  may  be  more  correct  in  terms  of  phys¬ 
ical  principles. 

Data  sources 

Most  of  the  data  used  In  1982  were  from  obscure  or  unpublished 
reports,  so  all  relevant  data  were  tabulated  in  an  appendix  to  the  report 
(Melior  1982).  Data  added  to  the  file  since  then  are  from  similarly 
obscure  sources,  and  they  are  given  here  in  Appendix  A.  Two  tables  give 
results  of  tests  on  arctic  sea  ice  that  were  made  in  1960  and  1961  by  the 
U.S.  Navy;  the  data  were  not  available  for  release  until  recently.  Another 
table  gives  some  results  from  tests  that  were  made  in  Korea  by  the  author 
subsequent  to  the  1982  report.  In  these  tests,  the  emphasis  was  on  pattern 
charges,  but  a  few  single-charge  shots  were  made.  Another  table  covers 
tests  made  by  the  author  at  a  site  in  the  U.S.  in  1984.  The  final  tables 
give  some  very  recent  results  from  Korea  (data  from  Table  A5  were  not  in¬ 
cluded  in  the  analysis). 


When  these  results  were  added  to  the  original  collection,  the  file 
that  was  analyzed  had  312  data  sets,  compared  with  291  in  1982  (315  sets 
are  now  available). 

Traditional  English  units  were  used  in  the  tabulation  and  analysis  of 
data  (App.  B),  since  these  units  still  dominate  in  North  American  explo¬ 
sions  technology.  Final  results  can  easily  be  converted  into  SI  units. 

Scaling  of  test  data  for  standard  curves 

To  remove  the  effect  of  charge  size,  Hopkinson  scaling,  or  cube  root 
scaling,  was  adopted  from  the  start.  All  linear  dimensions  of  the  problem 
are  divided  by  the  cube  root  of  charge  weight.  Although  this  gives  linear 
dimensions  with  units  of  ft/lb^/^  or  m/kg^^,  the  lengths  are  effect¬ 
ively  dimensionless,  since  the  cube  root  of  charge  weight  represents  the 
charge  radius  for  constant  explosive  density.  This  type  of  scaling  implies 
geometric  similarity. 

Cube  root  scaling  is  not  strictly  valid  for  underground  explosions, 
where  the  material  strength  increases  with  gravity  body  forces,  nor  for 
underwater  explosions,  where  the  static  head  is  the  water  depth  plus  a  10-m 
head  representing  atmospheric  pressure.  However,  in  neither  of  these  situ¬ 
ations  does  cube  root  scaling  distort  the  results  significantly  when  charge 
mass  ranges  from  less  than  1  kg  to  more  than  1  tonne.  To  improve  the  fit¬ 
ting  of  curves  to  cratering  data  or  water  plume  data  for  certain  ranges, 
empirical  exponents  smaller  than  one-third  have  sometimes  been  applied  to 
the  charge  weight  that  scales  either  depth,  or  all  linear  dimensions.  This 
is  inelegant  and  potentially  confusing;  it  is  best  avoided  if  possible. 

Alternative  normalizing  parameters  include:  (l)  charge  radius  for  an 
equivalent  spherical  charge  of  some  reference  explosive,  and  (2)  the  theo¬ 
retical  maximum  radius  of  the  underwater  gas  bubble  for  the  first  bubble 
pulse.  These  are  unsuitable  for  practical  work.  The  engineer  or  blaster 
knows,  or  can  determine,  the  weight  of  an  explosive  charge.  By  contrast, 
effective  charge  radius  and  maximum  bubble  radius  have  to  be  determined 
from  calculations  or  conversion  charts. 


Standard  data  analysis 

To  update  the  basic  design  curves,  the  original  analytical  procedure 
was  followed  (Mellor  1972,  1982).  Taking  scaled  charge  depth  and  scaled 
ice  thickness  as  independent  variables  Xj  and  X2,  the  scaled  crater  radius 
(dependent  variable  Y)  is  first  expressed  as  a  10-term  polynomial  in  Xj 


and  X2,  with  terras  and  cross-products  up  to  the  third  power.  Assuming  all 
error  to  be  in  determination  of  Y,  the  10  coefficients  of  the  regression 
equation  are  then  determined.  The  relative  significance  of  each  terra  is 
examined  by  stepwise  regression,  coefficients  being  deleted  successively 
from  t-test  values,  with  successive  checks  on  the  multiple  correlation  co¬ 
efficient  and  the  standard  error  of  estimate.  This  has  always  led  to  the 
dropping  of  the  same  two  terms,  bj  and  b4  of  the  10-term  polynomial  (coef¬ 
ficients  for  Xj  and  X[X2).  The  curves  are  then  drawn  from  this  final  8- 
term  regression  equation;  Y  is  plotted  against  Xj  with  X2  as  parameter,  and 
Y  is  plotted  against  X2  with  Xj  as  parameter  (Fig.  1). 

For  the  present  collection  of  data,  the  standard  error  of  estimate  is 
1.370  ft/lb^/3,  or  0.543  m/kg^/^  (a  deterioration  from  the  1982  value 
of  1.265  ft/lbl/3  or  0.502  m/kgl/3)  and  the  multiple  correlation  coef¬ 
ficient  r  =  0.708  (almost  the  same  as  the  1982  value  of  0.707).  The  8-term 
regression  equation  in  units  of  ft  and  lb  was: 

Y  =  5.222  +  3.208  X2  -  0.1273  X*  -  2.242  X2  +  0.006870  Xj 
+  0.02191  XjX2  -  0.07376  XxX2  +  0.3627  X2  3 


Effects  of  explosive  type  and  ice  type 

With  only  a  limited  amount  of  data,  it  is  not  feasible  to  consider  the 
variation  of  explosive  type,  or  of  ice  type,  in  the  regression  analysis. 

The  specific  energy  of  typical  explosives  varies  within  fairly  narrow 
limits,  and  there  is  appreciably  less  variation  in  the  cube  root  of  the 
charge  yield  that  is  used  for  scaling  linear  dimensions  (typically  <  ±9%). 
Only  a  small  proportion  of  the  total  energy  goes  into  shock  propagation, 
and  therefore  changes  in  that  proportion  relative  to  the  energy  for  gas 
expansion  are  not  of  major  significance.  In  1982  the  data  were  examined 
carefully  for  signs  of  systematic  variation  attributable  to  explosive  type, 
but  there  was  no  conclusive  evidence  of  such  an  effect.  Since  then,  it  has 
been  found  that  gas  blasting  devices  produce  much  the  same  ice-breaking 
effects  as  energetically  equivalent  charges  of  high  explosive  (Mellor 
1984a),  so  it  3eems  likely  that  variation  of  explosive  characteristics  is 
not  of  major  significance  when  the  degree  of  fragmentation  is  not  of  great 
concern. 

Data  for  underground  explosions  show  that  crater  dimensions  in  a  semi¬ 
infinite  solid  are  rather  insensitive  to  material  properties  for  a  wide 
range  of  rocks  and  soils,  and  this  suggests  that  minor  variation  of  ice 


idms  (ft/lb 


Scaled  Ice  Thickness  (ft/lb3) 

a.  Scaled  crater  radius  as  a  function  of  scaled  ice 
thickness,  with  scaled  charge  depth  as  parameter. 


b.  Scaled  crater  radius  as  a  function  of  scaled 
charge  depth,  with  scaled  ice  thickness  as  parameter 


Figure  1.  Design  curves  for  ice  blasting. 


properties  may  have  Little  effect  on  craters  in  ice-  In  a  previous  report 
(Mellor  1982),  it  was  suggested  that  the  thickness  of  a  snow  cover  on 
floating  ice  might  have  more  significance  than  variation  in  ice  properties. 

Interpretation  and  application  of  design  curves 

The  standard  curves  are  generated  in  traditional  English  units  for  use 
in  the  U.S.  Equivalent  SI  plots  can  be  produced  by  suitable  choice  of  in¬ 
tervals  for  the  parameter,  but  there  is  no  real  need  if  appropriate  conver¬ 
sion  factors  are  used  in  working  a  problem.  To  convert  ft/lb^-/^  to 
m/kgl/^,  multiply  by  0.3967;  to  convert  rn/kg^/-^  to  ft/lb^/^,  multiply 
by  2.5208.  For  ra-ft  and  kg-lb  conversion,  the  factors  are  3.281  and  2.205 
respect ive ly . 

Figure  1  shows  that  maximum  crater  diameter,  or  area,  is  achieved  with 
a  scaled  charge  depth  of  0  to  0-5  ft/lb^/^  (0  to  0.2  m/kg*-/^)  below  the 
base  of  the  ice  cover,  provided  that  the  scaled  ice  thickness  is  around  0-9 
ft/lb'-/^  (0.36  m/kg*^^).  In  other  words,  the  best  results  are  obtained 
with  the  charge  almost  in  contact  with  the  base  of  the  ice.  With  these 
optimum  conditions,  the  scaled  crater  radius  is  6-56  ft/lb*-/^  (2-6 
m/kg  The  predicted  scaled  crater  diameter  is  thus  just  over  13 

fr/lb1/^  (5.2  m/kg1/^),  or  about  15  times  the  ice  thickness  (Fig.  2). 

For  scaled  charge  depths  greater  than  0.5  ft/lb1/^  (Q.2  m  /k g1/3) > 
scaled  crater  radius  decreases  fnr  all  values  of  scaled  ice  thickness 
(Fig.  1).  At  these  greater  charge  depths,  the  optimum  value  of  scaled  ice 
thickness  appears  to  vary  somewhat- 

When  using  the  curves,  the  usual  input  Is  a  measured  or  estimated 
value  of  the  ice  thickness  t.  If  the  intent  is  to  get  the  biggest  crater 
for  the  least  amount  of  explosive,  an  optimum  charge  weight  WQpt  can  then 
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Figure  2-  Simple  guidelines  for  optimum  charge 
design  (single  crater). 


be  estimated  from  the  optimum  ice  thickness,  which  is  approximately  0.9 
ft/lb1/3  (or  0.36  m/kg1/3): 

t/Wl/3  -  0-9  ft/lb1/3 

opt 

or 

W  -  (t/0.9) 3  lb  *  1.4t 3  lb 

opt  v 

when  t  is  in  ft.  If  t  is  in  m,  then 
W  -  (t/O-36)  3  «  21t  3  kg. 

The  optimum  charge  depth  dc/wl/3  0  to  0.5  ft/lb1/3  (0-0.2 

m/kg1/3) ,  so  that 

d  «  0  -  0 . 6t . 
c 

The  easiest  rule-of-thumb  for  predicting  the  diameter  D  of  an  optimum 
crater  Is  to  take 
D  -  1 5t  • 

With  ice  thickness  given,  the  curves  can  provide  predictions  of  crater 
size  for  non-optimum  charge  weight  and/or  charge  depth,  or  they  can  be  used 
to  determine  the  size  of  charge  needed  to  produce  a  crater  of  specified 
diameter  from  a  given  charge  depth.  Alternatively,  the  optimum  or  maximum 
thickness  of  ic.e  can  be  calculated  when  charge  size  and  charge  depth  are 
specified-  These  various  applications  have  been  illustrated  by  means  of 
worked  numerical  examples  (Me.l1.or  1982). 


Multiple  charges  for  ice  breaking 

A  row  of  under-tce  charges  can  be  used  to  break  a  channel,  and  a 
multi-row  array  can  be  used  to  fracture  ice  over  a  broad  area-  To  design  a 
row,  or  an  array,  of  ice- breaking  charges,  it  is  convenient  to  first  calcu¬ 
late  the  crater  radius  Rc  ^  for  a  single  charge  and  then  to  relate  Rc^ 
to  the  spacing  s  for  charges  in  a  row  or  in  a  square  network. 

In  a  single  row  of  charges  where  s/Rcj  is  greater  than  about  2.5,  or 
in  an  array  where  s/Rr]  >  4,  each  charge  produce®  a  separate  eruption  and 
a  separate  crater  (see  Table  l  for  references).  When  s/Rc^  becomes  some¬ 
what  smaller  r ban  these  values,  each  charge  still  gives  a  separate  eruption 
(fig.  3b- g),  hot  as  the  ice  debris  and  the  water  columns  fall  back  again, 
the  i ee  between  the  individual  craters  is  subjected  to  impact,  base  surge, 
and  violent  wave  action  CFtg.  3  h-m).  This  can  fracture  the  ice  between 
th>>  craters  with  varying  degrees  of  fragmentation,  depending  on  s/Rcj. 


Figure  3  (coat'd).  An  array  of  multiple  charges  breaking  through  thi 
ice  (3-lb  charges  on  27-ft  centers).  Charge  depth  0.45  bubble  radii. 
Photos  in  this  sequence  were  taken  at  intervals  of  approximately  0.2 
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Figure  3  (coat'd).  An  array  of  multiple  charges  breaking  through  thin 
ice  (5-lb  charges  on  27-ft  centers).  Charge  depth  0.45  bubble  radii. 
Photos  in  this  sequence  were  taken  at  intervals  of  approximately  0.2  s 


Table  1.  Effects  of  charges  in  a  single  row  when  each  charge  is  close  to 
optimum  for  the  prevailing  ice  thickness,  and  charge  depth  is  about  It  to 
1.5t.  Based  on  tests  by  Kurtz  et  al*  (1966),  Fonstad  and  Gerard  (1985)  and 
unpublished  work  by  the  author. 

Approximate  width  of 
Charge  spacing  fractured  channel,  b 


(s/Rcl) 

(b/Rcl) 

(b/s) 

Fragmentation 

>2.5 

Individual 

craters 

As  for  single  crater. 

2 

3.8 

1.9 

Poor  at  edges  of  chan¬ 
nel  and  in  cusps  be¬ 
tween  craters. 

1.5 

3.4 

2.3 

Good  in  mid-channel, 
less  along  edges. 

1 .0 

3.4 

3.4 

Good. 

Under  these  conditions,  the  author  has  chosen  to  accept  1  ft/lb1/3  (0.4 
m/kg1/3)  as  optimum  charge  depth  (instead  of  0-0.5  ft/lb1/3)  in  order 
to  increase  the  base  diameter  of  the  water  plume.  There  are  no  systematic 
data  to  support  this  judgment,  which  derives  from  findings  for  explosions 
beneath  ice-free  surfaces  (see  MelLor  1986). 

Drawing  upon  limited  field  data,  probable  effects  of  row  charges  are 
indicated  in  Table  1.  The  charges  interact  to  produce  a  continuous  broken 
channel  when  s/Rc^  <  2. 5-  Fragmentation  improves  and  becomes  more  uni¬ 
form  as  s/Rci  decreases,  until  the  row  of  individual  charges  approximates 
the  effect  of  a  continuous  linear  charge,  for  which  the  blast  effects 
spread  cylindrical ly  rather  than  spherically. 

A  multi-row  array  produces  more  interference  between  craters,  and  the 
charge  spacing  can  be  wider  than  it  is  in  a  single  row.  Hie  logical  charge 
pattern  is  a  square  net  of  mesh  size  s.  A  "diamond,"  or  "5-spot,"  pattern 
is  obtained  by  rotating  the  reference  axes  of  a  square  mesh  through  45°. 
Limited  personal  experience  with  charge  arrays  indicates  that  Individual 
craters  are  formed  with  s/Rc^  >  4,  continuous  fracture  occurs  with 
s/Rcj  =  3,  adequate  fragmentation*  is  obtained  with  s/Rc^  =  2-7, 
spacing  is  somewhat  closer  than  necessary  with  s/Rc^  =  2-3,  and  spacing 
is  wasteful ly  close  with  s/Rcj  =  2.  Fonstad  and  Gerard  (1985)  deduced  an 

*More  or  less  equant  fragments,  as  distinct  from  slab  fragments. 


upper  limit  for  useful  interference  ns  s/Rc  j  =  3 ■ S ,  which  agrees  reason¬ 
able  well  wi*l>  <  / 1  i'.u*  of  1.8  suggested  earlier  (Me)lor  198?). 

Charge  design  for  thin  ice 

It  the  ice  's  thin  Ml  ft,  ' .  1  m)  ,  the  optimum  size  of  a  single 
charge  is  small,  and  the  absolute  diameter  of  the  resulting  crater  is  also 
small  ( <  l  a  ft,  <4  .  b  m)  •  When  tlie  objective  is  to  break  a  wide  channel  or  a 
broad  area  in  thin  ice  wit'1  multiple  charges,  the  required  number  of  opti¬ 
mum  charges  may  be  Large.  Under  these  circumstances,  the  effort  involved 
in  drilling  and  charge  placement  can  become  excessive,  and  it  may  he  more 
efficient  in  practical  terms  to  use  charges  that  are  bigger  than  optimum, 
and  to  set  them  further  apart. 

Having  made  the  calculations  for  opt imem  charge  size  to  get  a  fee)  for 
th'ogs,  guess  can  he  mads  at  the  total  amount  of  explosive  that  seems  af 
fo'dtble  for  the  job,  and  a  new  charge  size  on  be  estimated  (e-g.  5  lb  in¬ 
stead  of  1  lb,  or  2.5  kg  1  *ist  end  of  n.  s  kg'  The  scaled  ice  thickness  for 
the  new  charge  weight  is  then  calculated.  and  the  curves  of  Figure  1  are 
used  to  obtain  the  corresponding  value  of  Pr j  (perhaps  taking  a  scaled 
charge  depth  of  abet  1  ft/lh'^  or  0  ■  h  m/kg 1/1  -  ,  as  suggested  ear¬ 
lier).  'icjng  the  chosen  value  of  *7  P.c\ >  say  2.7,  the  new  value  of  charge 
spacing  s  is  obtained- 

Specific  energy,  or  powder  factor,  for  ice  fragmentation 

The  specific  energy  for  f ragmen  fat ion  of  a  solid,  Es ,  1 s  the  energy 
consumed  to  break  unit  volume.  The  energy  in  a  charge  of  weight  W  is  kW, 
where  k  ica  /'h<a  rac  t  e  r  i  s  t  i  c  specific  energy  content  for  the  particular  ex¬ 
plosive.  k  ran  be  taken  either  as  the  heat  of  explosion,  typically  around 
U  -  2  to  ,  q  k.1 7  a  for  high  explosives,  or  as  the  gas  expansion  energy,  which 
is  roiighlv  ?.U%  of  the  heat  of  explosion.  The  volume  of  Ice  broken  hv  a 
single  charge  is  R,,  t ,  where  Rc  is  crater  radius  and  t  is  ice  thick¬ 
ness.  This,  for  explooi-.e  icebreaking, 
kW  k/  w 

The  weight  of  ° x p  1  o s  I  .  ■  p  jer  uni  t  volume  of  fragmented  material  Is 
known  as  the  powder  factor  F  or  as  the  specific  charge-  For  ice  Mast¬ 
ing  with  a  single  r  h  a  r  g  e : 


1  7  ,  1/1 
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Figure  4.  Specific  charge,  or  powder  factor,  as  a 
function  of  scaled  ice  thickness  when  single  charges 
are  at  optimum  depth. 


Figure  4  gives  the  specific  charge  (powder  factor)  as  a  function  of 
scaled  ice  thickness  for  single  charges  at  optimum  depth,  i.e.  immediately 
below  the  bottom  of  the  ice-  The  minimum  value  of  specific  charge  is 
0-0051  lb/ft1,  or  0.082  kg/nr,  and  it  is  obtained  when  the  scaled  ice 
thickness  is  about  2  ft/lb*/-*,  or  0.8  m/kg*/^.  In  other  words,  the 
maximum  specific,  volume  of  fragmentation  occurs  with  t/W^V-^  *  2 
ft/ib*/^  (0-8/m  kg^/^),  whereas  maximum  specific  area  of  fragmentation 
occurs  wi  tb  t/wl/3  a  0.9  ft/lb^^  (0-36  m/kg^/-^)- 

For  comparison,  the  powder  factor  for  optimum  crater  blasts  in  semi¬ 
infinite  ice  is  about  0.01  to  0.02  lb/ft^,  or  0.004  to  0*008  kg/m^. 

Values  of  specific  energy  Es  for  explosive  icebreaking  by  single 
charges  can  be  obtained  from  Figure  4  by  multiplying  values  of  specific 
charge  by  either  the  heat  of  explosion  or  the  energy  of  gas  expansion  for 
explosive.  The  units  of  this  specific  energy  are  energy  per  unit  volume, 
which  reduces  to  the  dimensions  of  a  stress.  If  we  take  the  heat  of  explo 

sion  as  4.6  kJ/g  (l.l  kcal/g),  which  is  about  right  for  TNT,  the  minimum 

3 

overall  specific  energy  for  explosive  icebreaking  is  about  0.38  MJ/m  , 
which  is  equivalent  to  0-38  MPa  (55  lbf/in.^).  Tills  is  comparable  to  the 


specific  energy  of  a  fairly  large  icebreaker  (10,000-50,000  hp)  when  Es 
is  based  on  shaft  horsepower.  The  "process  specific  energy”  for  a  ship, 
based  on  propeller  thrust  and  speed,  is  lower,  but  this  should  be  compared 
with  an  explosive  specific  energy  that  is  based  on  the  work  of  gas  expan¬ 
sion  (about  20%  of  the  value  given  above). 

The  specific  energy  for  maximum  crater  area  is  obtained  by  taking 
rc/w1/3  =  6.56  ft/lb1/3  and  t/W1/3  =  0.9  ft/lb1/3.  These  give 
Fp  =  0.0082  lb/ft3  (0.13  kg/m3)  and,  with  k  =  4.6  kJ/g,  Es  =  0.6 
MJ/m3  (0.6  MPa,  88  lbf/in.2).  If  specific  energy  is  based  on  the  work  of 
gas  expansion,  Eg  »  0.12  MJ/m3.  These  values  are  a  bit  higher  than  cor¬ 
responding  values  for  maximum  specific  volume  of  broken  ice,  but  are  still 
comparable  to  the  specific  energy  for  an  icebreaking  ship. 

Other  analytical  schemes 

In  practical  ice  blasting,  charge  depth  has  usually  been  referred  to 
the  base  of  the  ice  sheet,  and  test  data  have  been  reported  accordingly. 
However,  in  underwater  explosion  technology,  the  water  surface  is  the 
obvious  datum  for  charge  depth,  and  it  can  be  argued  that  it  is  also  a 
logical  datum  for  charge  depth  in  ice  blasting.  To  convert  traditional 
charge  depth  to  charge  depth  below  water  level,  we  simply  add  the  ice 
thickness  multiplied  by  the  specific  gravity  of  the  ice.  Having  made  this 
conversion,  the  regression  analysis  can  be  repeated  to  give  a  prediction  of 
scaled  crater  radius  as  a  function  of  scaled  ice  thickness  and  scaled 
charge  depth  below  water  level. 

Figure  5  gives  the  prediction  curves  when  water  level  is  the  reference 
for  charge  depth.  The  curves  are  broadly  similar  to  the  standard  curves, 
but  they  necessarily  cut  off  for  the  portions  where  depth  is  less  than  the 
ice  thickness.  The  predictions  are  essentially  the  same  as  those  given  by 
the  standard  curves.  However,  the  empirical  correlation  is  not  as  good  as 
it  is  for  the  standard  curves.  The  standard  error  of  Y  is  1.408 
ft/lb^-/3,  compared  with  1.370  ft/lb^^3.  The  multiple  correlation  coef¬ 
ficient  r  is  0.684,  compared  with  0.708  for  the  standard  curves. 

In  Figure  5b  the  ice  cratering  behavior  is  compared  with  the  types  of 
surface  eruptions  which  would  occur  at  an  ice-free  surface.  The  maximum 
scaled  crater  radius  is  achieved  when  charge  depth  is  in  the  shallow  end  of 
the  range  which  produces  columnar  eruptions  from  a  free  surface. 


Chorqe  Deo’fi  Beio*  Wo ter  Level  (bubble  ro<JnJ 

figure  6.  Prediction  curves  derived  from  regres 
sion  analysis  when  linear  dimensions  are  scaled 
with  respect  to  maximum  bubble  radius,  and  water 
level  is  the  depth  datum. 


In  underwater  explosion  technology,  linear  dimensions  are  sometimes 
scaled  with  respect  to  the  theoretical  maximum  radius  of  the  gas  bubble  for 
the  first  pulse,  R^m’  This  procedure  is  particularly  useful  in  dealing 
with  problems  that  involve  surface  effects,  such  as  eruption  of  water¬ 
spouts.  Using  standard  relations  which  give  Rjjm  as  a  function  of  charge 
weight  and  depth  below  the  water  surface  (see  Mellor  1986),  the  data  for 
ice  blasting  can  be  converted  into  dimensionless  form,  such  that  all  linear 
dimensions  are  normalized  with  respect  to  R^*  The  data  were  scaled  in 
this  way,  using  water  level  as  the  depth  datum,  and  the  regression  analysis 
was  repeated- 

Results  of  the  regression  analysis  for  the  bubble-scaled  data  are 
shown  in  Figure  6.  The  predictions  are  more  or  less  the  same  as  those  of 
the  standard  curves,  but  the  empirical  correlation  is  not  as  good.  The 
multiple  correlation  coefficient  r  is  0.649,  compared  with  0.708  for  the 
standard  curves  and  0-684  for  cube-root  data  referred  to  water  level.  In 
Figure  6b  the  depth  ranges  for  various  types  of  free-surface  eruptions  are 
shown.  As  before,  it  can  be  seen  that  the  scaled  crater  radius  is  a  maxi¬ 
mum  when  charge  depth  is  in  the  range  that  produces  columnar  waterspouts  in 
ice-free  conditions. 

For  practical  use,  there  is  no  doubt  that  the  standard  curves  are  more 
convenient,  and  slightly  more  reliable,  than  the  alternative  forms  de¬ 
scribed  here.  However,  the  bubble-scaled  curves  provide  some  insight  into 
the  relation  between  ice  cratering  and  underwater  gas  bubbles.  The  maximum 
scaled  crater  radius  that  can  be  achieved  is  about  1.63  Rjjm,  and  it  is 
obtained  when  the  charge  depth  and  the  ice  thickness  are  only  about  0.3 
R^m.  This  implies  that  water  driven  by  the  expanding  gas  bursts  through 
the  ice  during  the  first  bubble  expansion  (Fig.  7). 

Gas  blasting 

Floating  ice  can  be  broken  by  a  rapid  gas  discharge  beneath  the  ice. 
Tests  have  been  made  with  a  variety  of  systems,  all  of  which  discharge  at 
relatively  low  pressure  without  propagating  a  true  shock  (Mellor  1984a). 
Experimental  equipment  has  included:  (a)  carbon  dioxide  shells  (discharge 
pressure  70-90  MPa),  (b)  an  airblasting  system  (69-83  MPa),  (c)  an  air  gun 
(17  MPa),  and  (d)  a  fuel/air  combustion  system  (2. 3-4. 3  MPa). 

Gas  blasting  results  cannot  be  compared  directly  with  explosives 
data.  To  draw  any  kind  of  conclusions  from  the  limited  gas  blasting  data, 
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Figure  7.  Eruption  through  very  thin  ice 
below  water  level  is  0.49  bubble  radii, 
'■motion  is  a  column  with  a  center  jet. 


cl  ime  ns  ions  have  to  be  scaled  with  respect  to  the  cube  root  of  the  discharge 
energy,  which  is  taken  as  the  energy  for  adiabatic  gas  expansion  to  atmo¬ 
spheric  pressure.  For  comparison  with  explosives,  it  is  not  clear  whether 
the  energy  of  the  explosive  should  he  taken  ;»s  the  heat  of  explosion  or  as 
the  energy  of  gas  expansion.  If  a  scaled  dimension  for  gas  blasting  has 
units  of  it  can  he  converted  to  units  of  m/kg^/^  hy  multiply¬ 

ing  by  the  cube  root  of  an  explosive  specific  energy  factor  that  has  units 
of  MJ/kg.  For  TNT,  the  heat  of  explosion  is  approximately  4.56  MJ/kg  and 
the  energy  of  gas  expansion  is  0.87  MJ/kg;  the  corresponding  cube  root  fac¬ 
tors  are  1.66  and  0.955  (MJ/kg)^/-^. 

For  ice-breaking  gas  discharges  of  magnitude  L  MJ ,  optimum  scaled  ice 
thickness  for  maximum  crater  diameter  is  about  0.4  optimum 

charge  depth  is  in  the  range  0  to  0.6  m/MJ '■/■*,  and  the  resulting  maximum 
crater  radius  is  about  2.9  (Mellor  1984a).  If  the  gas  expansion 

conversion  factor  tor  high  explosives  is  applied  to  the  optimum  crater 
results  for  explosives  (Fig.  1),  the  corresponding  values  are:  optimum  ice 
thickness  0.38  optimum  charge  depth  0  to  0.2.1  an(j 

maximum  crater  radius  2.7  Considering  all  the  uncertainties 

that  are  involved,  this  is  remarkably  close  agreement  between  gas  blasting 
and  explosives  blasting.  Furthermore,  both  gas  blasting  and  explosives 
give  a  probable  maximum  value  for  the  scaled  crater  diameter  of  approxi¬ 
mately  15  times  the  optimum  ice  thickness. 

The  specific  energy  Es  for  maximum  values  of  s.  aled  crater  diameter 
is  about  0.1  MJ/m^.  This  is  the  same  as  the  corresponding  specific  energy 
for  conventional  explosives  when  Es  is  based  on  the  work  of  gas  expansion 
rather  than  upon  the  heat  of  explosion  ( Es  ~  0.12  MJ/m  )• 

This  suggests  that,  for  optimum  cratering,  the  crater  size  is  con¬ 
trolled  by  gas  expansion  and  water  displacement,  not  by  shock  propagation. 
Since  most  of  the  energy  of  a  high  explosive,  goes  into  gas  expansion,  and 
since  the  energy  per  unit  mass  of  typical  explosives  does  not  ''ary  between 
very  wide  limits,  It  also  seems  reasonable  to  expect  that  optimum  '’rater 
dimensions  will  be  fairly  Insensitive  to  explosive  type. 

if  crater  size  is  influenced  mainly  by  gas  expansion  and  water  dis¬ 
placement,  then  variations  in  ice  type  might  not  have  much  effect*  In  par¬ 
ticular,  an  Ice  cover  formed  by  accumulation  of  ice  fragments  might  respond 
much  like  an  intact  ice  plate. 


Charges  on,  or  within,  the  ice  sheet 

A  small  charge  placed  on  top  of  a  thick.  Ice  sheet,  or  in  its  upper 
half,  has  much  the  same  effect  as  a  cratering  charge  in  semi-infinite  ice 
(Fig-  8).  Figure  9  summarizes  the  data  for  scaled  dimensions  of  the  true 
crater  in  deep  ice-  The  true  crater  is  the  more  or  less  conical  region  of 
fractured  ice,  as  distinct  from  the  apparent  crater,  which  is  the  visible 
open  hole.  The  maximum  scaled  radius  of  the  true  crater,  measured  at  the 
ice/air  surface,  is  in  the  range  3  to  5  ft/lb*/-*  (1-2  to  2-0  m/kg^/^). 

This  maximum  is  achieved  with  a  scaled  charge  depth  of  3  to  4  ft/lb^^ 

(1.2  to  1-6  m/kg1/^). 

Corresponding  data  for  the  apparent  crater  are  given  in  Figure  10. 
Dimensions  of  the  apparent  crater  are  smaller  than  those  of  the  true  crater 
because  some  of  the  fractured  ice  is  not  grossly  displaced,  and  some  ejec¬ 
ted  fragments  fall  back  into  the  crater. 


(m/kq  )  fft/lb 


Scaiert  Chorqe  Depth 

b.  Scaled  crater  depth  as  a  function  of  scaled 
charge  depth. 

Figure  9.  Summitry  at  data  tor  the  scaled  dimen 
sions  of  the  true  crater  in  massive  ice.  (Data 
from  Livingston  I960.) 
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a.  Scaled  crater  radius  as  a  function  of 
scaled  charge  depth. 


b.  Scaled  crater  depth  as  a  function  of 
scaled  charge  depth. 


Figure  10.  Summary  of  scaled  dimensions 
for  apparent  craters  in  massive  ice. 
(Data  from  Livingston  I960.) 


*1  Small  cratering  charge  to 
control  hole  diameter 
* 2  Larger  charge  venting 

Charges  tired  in  sequence  with  delays  through  crater  ol  «  1 

in  order  to  produce  narrow  hole 

and  to  clear  deOns  from  hole  '  3  Cfta,9e  ''en,mS  uP*a'd  a"d 

also  breaking  to  base  of  ice 

*4  Relatively  large  clearing 
charge  below  ice 


Figure  11.  Use  of  delay  deck  charges  to  break  a 
narrow  shaft  through  very  thick  ice. 

A  charge  set  at,  or  below,  mid-depth  can  break  out  to  both  the  ice/air 
and  the  ice/water  interfaces,  provided  that  t/2W^/3  is  less  than  about  4 
ft/lb1/3  (1.6  m/kg1/3).  However,  the  maximum  scaled  radius  of  the 
upper  or  lower  surface  fracture  zone  is  likely  to  drop  below  the  3-5 
ft/lb1/3  (i. 2-2.0  m/kg1/3)  value  for  an  optimum  cratering  charge  in 
very  thick  ice. 

To  cut  a  relatively  narrow  vertical  shaft  through  thick  ice,  delayed 
deck  charges  can  be  employed-  Several  charges  are  placed  at  different 
depths,  and  they  are  fired  in  rapid  sequence,  from  top  to  bottom  with  ap¬ 
propriate  delays  (Fig.  11).  The  top  charge  is  designed  to  produce  the 
required  hole  diameter,  in  accordance  with  the  cratering  data  of  Figure  9. 
Successive  charges  vent  through  this  crater,  and  a  final  charge  in  the 
water  beneath  the  ice  flushes  out  ice  fragments  (which  unfortunately  tend 
to  wash  back  into  th)  hole  as  water  drains  from  the  ice  surface). 

When  a  small  unconfined  charge  is  laid  on  top  of  thick  ice  it  makes  a 
superficial  crater,  with  the  dimensions  given  by  Figures  9  and  10  for  zero 
charge  depth.  The  scaled  radius  of  the  true  crater  is  2-3  ft/lb1/3  (0.8 
-  1.2  m/kg1/3)  at  the  surface,  but  the  depth  is  small,  say  0.6  to  1.0 
ft/lb1/3  (0.24  to  0.4  m/kg1/3).  In  order  for  a  surface  charge  to  break 
through  to  the  bottom  of  the  ice  sheet  and  to  form  a  penetration  crater,  as 
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Figure  12.  Effect  of  a  charge  lying  on  the  upper 
surface  of  thin  ice. 


Figure  13.  Effects  of  surface 
charges  compared  with  the  effects 
of  under-ice  charges.  [Data  for 
surface  charges  from  Barash  (1966) 
and  Fonstad  et  al.  (1982).] 


in  Figure  12,  the  scaled  ice  thickness  has  to  be  less  than  2  ft/lb^^ 
(0.8  m/kgl/3).  As  the  scaled  ice  thickness  decreases,  the  radius  of  the 
crater  increases,  but  it  is  always  much  smaller  than  the  radius  of  the 
crater  that  would  be  produced  by  the  same  charge  set  in  water  under  the 
base  of  the  ice  sheet  (Fig.  13).  A  surface  charge  is  likely  to  produce 
more  breakage  if  it  is  covered  by  sandbags  or  suchlike  ("mudcapping"  or 
"dobying”) . 
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A  charge  la  air  above  the  Ice  surface  is  even  less  effective  than  a 
surface  charge-  No  systematic  test  results  have  been  reported,  but  at  a 
scaled  ice  thickness  of  1  ft/lb^/^  (0.4  m/kg^-/^  ) ,  no  crater  is  formed 
when  the  charge  height  is  0.5  ft/lb^/^  (0.2  m/kg^/^)  or  more. 

Pene t ration  by  shaped  charges 

A  shaped  charge,  or  lined  cavity  charge,  forms  a  high  velocity  jet  and 
incorporates  into  it  dense  material  from  a  metal  or  glass  cavity  liner. 
Typical  charges,  which  are  rotationally  symmetrical,  produce  a  narrow  hole 
in  the  target  material.  There  are  also  linear  shaped  charges,  which  pro¬ 
duce  a  slot  in  the  target. 

Some  guidelines  for  ice  penetration  by  typical  shaped  charges  with 
conical  liners  have  been  developed  (Mellor  1984b),  although  test  results 
are  few.  Depending  on  the  cone  angle  of  the  charge,  penetration  depth  in 
thick  ice  is  likely  to  be  in  the  range  12  to  16  cone  diameters,  with  mean 
hole  diameter  in  the  range  1/3  to  2/3  of  the  cone  diameter.  Deep  penetra¬ 
tion  (16  cone  diameters)  and  a  slender  hole  (1/3  cone  diameter)  are  likely 
to  be  achieved  with  narrow  cone  angles,  say  around  45°.  Smaller  penetra¬ 
tion  (12  cone  diameters)  and  a  wider  hole  are  likely  with  wide-angle  cones 
(>  60°).  Optimum  standoff  distance  is  likely  to  be  in  the  range  2  to  4 
cone  diameters,  depending  somewhat  on  cone  angle.  All  of  these  values  are 
for  a  shaped  charge  detonating  in  air  above  the  ice  (Fig.  14). 
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P  greatest  with  cone  angle  42° -45°.  S  *  2  Dc -4  Dc 
Dh  greatest  with  cone  angle  60°  -**80° 

P  *  16  Dc  with  narrow  angle  cone 


For  complete  penetration.  D  >  t/16 


A  shaped  charge  fired  underwater  needs  a  standoff  cavity  where  the  jet 
can  develop  before  it  strikes  the  target  material.  This  can  be  a  canister 
filled  with  air  or  with  low-density  foamed  plastic.  Penetration  by  the 
explosive  jet  should  then  be  similar  to  the  penetration  from  air,  but  con¬ 
tainment  of  the  charge  by  water  could  lead  to  greater  local  cratering  in 
the  ice  surface. 

Conclusion 

Empirical  prediction  curves  for  explosive  ice-breaking  have  been 
developed  in  various  forms  from  test  data,  using  regression  analysis.  Hie 
"standard"  form  given  in  Figure  1  is  convenient  for  practical  use,  and  it 
appears  to  represent  the  data  slightly  better  than  alternative  forms,  which 
may  have  somewhat  closer  relation  to  the  relevant  physics  (Fig.  5  and  6) • 
There  is  still  likely  to  be  considerable  error  when  the  result  of  a 
particular  blast  is  compared  with  the  statistically  based  prediction.  How¬ 
ever,  the  present  guidelines  are  an  improvement  over  some  earlier  rules  for 
ice  blasting,  which  ignored  ice  thickness  as  a  variable  and  gave  optimum 
charge  depth  as  an  absolute  value,  irrespective  of  charge  size. 

The  main  requirement  for  improving  prediction  capability  is  more  data 
from  well-designed  field  tests,  including  tests  with  multiple  charges,  sur¬ 
face  charges,  charges  above  the  ice,  linear  charges  above  and  below  the 
ice,  and  shaped  charges  above  and  below  the  ice.  Test  shots  should  include 
conditions  that  are  far  from  optimum,  and  crater  dimensions  should  be 
defined  consistently. 
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APPENDIX  A:  TEST  DATA  ADDED  TO  THE  FILE  SINCE  PREVIOUS  REPORT  (MELLOR  1982) 


Table  A1 •  Project  Ice  Skate,  Bering  Sea,  February  1960.  U.S.  Naval  Ordnance 
Laboratory  (now  U.S.  Naval  Surface  Weapons  Center)*.  Explosive:  HBX-3. 


Charge  weight 
(lb) 

Average  ice 
thickness 
(ft) 

Depth  of  charge 
below  base  of  ice 
(ft) 

Average  crater 
diameter 
(ft) 

630 

4.6 

30 

100 

630 

2.3 

20 

no 

630 

2.6 

45 

90 

♦Leslie  and  Nelson  (1961).  Data  and  report  provided  by  courtesy  of  Robert  M 
Barash,  Explosion  Damage  Branch,  USNSWC. 


Table  A2 .  Polar  Ice  Experiment,  Beaufort  Sea,  Summer  1961.  U.S.  Naval  Ord¬ 
nance  Laboratory  (now  U.S.  Naval  Surface  Weapons  Center)*.  Explosive:  HBX-3. 


Charge  weight 

...Jib). _ 


Average  ice 
thickness 


(ft). 


Depth  of  charge 
below  base  of  ice 
.  (ft) 


Average  crater 
diameter 
_  (ft)_ _ 


60 

35 

20 


10  3.91  100 

14  3.27  40 

9.3  2.71  0 


♦Leslie  and  Nelson  (1961).  Data  and  report  provided  by  courtesy  of  Robert  M 
Barash,  Explosion  Damage  Branch,  USNSWC. 


Table  A3.  Tests  by  M.  Mellor,  Imjin  River,  Korea,  February  1983.  Water  depth: 
20  ft  (6  ra) ,  explosive:  C-4. 


Charge  weight 


Charge  depth 


Ice  thickness 


Crater  diameter 


Table  A4  .  Tests  by  M.  MelLor  at  Bolt  Technology  test  site,  January  1984.  Ex¬ 
plosive:  AN'  +  NM  binary 


Charge  weight 


Charge  depth 
below  base  of  ice 


Ice  thickness 


Crater  diameter 


Table  A5-  Tests  made  on  Imjin  River,  Korea,  by  E  Company,  2nd  Engineer 
Battalion,  U.S.  Army  (these  results  not  included  in  analysis) - 


Explosive 

C4 

Military 
dynami te 


Charge 

weight 

(lb) 

Charge  depth 
below  base  of  Ice 
(ft) 

Ice 

thickness 

(in) 

Crater 
diamete  r 
(ft) 

1.25 

0.33 

9 

18 

1.0 

0.33 

11.25 

8 

2.5 

0.33 

11 

20 

Table  A6  •  Pattern  shot  on  Imjin  River,  Korea,  by  E  Company,  2nd  Engineer  Bat 
talion,  U.S.  Army  (data  recorded  by  D.  Calkins,  CRRE1,,  January  1986.) 


Explosive 

C4 


Water  depth 

6  -  8~Tt 


Weight  of 
single  charge 

1.25  lb  (2.5  lb  in  center  row) 

Ice  thickness 
(average) 

0.89  ft 


Charge  depth 
0.33  ft  " 

Pattern 
Square  mesh, 

5  rows  18  ft  apart, 
16  charges  per  row 
at  18-ft  centers 


“cl 

6.52  ft/ lb1-/ 3 

s/Rci 

2.1b  ' 


Fragmen ta  t ion 
Equant  fragments  in 
central  crater  zones, 
some  small  slabs  near 
midpoints  of  charge  net. 
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APPENDIX  B:  CURRENT  DATA  FILE 


Test  data  are  tabulated  In  the  following  columns: 

-  charge  depth  below  base  of  ice  sheet  (ft/lb^-/^) 

X2  -  ice  thickness  (in./lb1/^) 

X3  -  ice  thickness  (ft/lb^/-*) 

X4  -  water  level  to  base  of  ice  (ft/lb^-^) 

X5  -  charge  depth  below  water  level  (ft/lb*/^) 

Xg  -  charge  weight  (lb) 

X7  -  theoretical  maximum  bubble  radius  during  first  pulse  (ft/lb^/-*) 
Xg  -  charge  depth  below  water  level  scaled  with  respect  to  X7 
Xg  -  ice  thickness  scaled  with  respect  to  X7 
Y1  -  crater  radius  (ft/lb1-/^) 

Y2  -  crater  radius  scaled  with  respect  to  maximum  bubble  radius,  X7 


The  actual  measured  values  are  scaled  and  entered  as  Yj,  Xj  and  X2/X3- 
The  charge  weight  is  listed  as  Xg.  All  other  quantities  are  generated  from 
these  values. 
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